#ifndef DYN_QUAD_QGAUSS4_H
#define DYN_QUAD_QGAUSS4_H

#include "dyn_quad_qGauss.h"

/*! \author Joey Dumont <joey.dumont@gmail.com>
 *
 *  \date 2012-10-07
 *
 *  \brief Gaussian quadrature with \f$N=4\f$.
 *
 * We declare the weights and abscissas associated
 * with \f$N=4\f$ Gaussian quadrature. The
 * values are taken from Abramowitz and Stegun.
 */

class qGauss4 : public qGauss
{
public:
    /*! Constructor defines functions and
     * boundaries.
     */
    qGauss4(Functor& _func, double _a, double _b);

    /*! Populate the vectors. */
    void populateVectors();

protected:
    /*! @name Abscissas */
    //@{
    double x1 = 0.339981043584856;
    double x2 = 0.861136311594053;
    //@}

    /*! @name Weights */
    //@{
    double w1 = 0.652145154862546;
    double w2 = 0.347854845137454;
    //@}
};

#endif // DYN_QUAD_QGAUSS4_H
